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SUMMARY 

A  mathematical  model  of  the  Sea  King  Mk.  50  anti-submarine  warfare  ( AS IV) 
helicopter  and  its  sonar  system  is  presented.  The  model  represents  both  performance 
and  dynamic  flight  behaviour  over  a  range  of  conditions  and  incorporates  the  aerodynamics 
and  kinematics  of  the  helicopter,  the  control  systems,  pilot  inputs,  the  cable  and  sonar 
dynamics,  plus  wind  and  sea  state  data.  The  aerodynamics  and  kinematics  model  is  a 
three-dimensional  representation  covering  the  operating  flight  envelope,  where  rotor 
aerodynamics  are  based  on  blade-eleir  •  /  actuator-disc  theory,  and  the  control 
systems  contain  models  of  both  the  fix , .  o/s  and  automatic  flight  control  system. 

The  sonar  cable  and  transducer  mode!  .\>  ,u,  '  by  a  number  of  attached,  rigid  links, 

where  consideration  of  the  forces  acting  on  cadi  link  enable  the  three-dimensional  shape 
and  motion  of  the  complete  cable  and  transducer  to  be  predicted.  A  description  of  the 
computer  program  for  the  matlnmatical  mode!  is  given  and  sample  results  are  included. 
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ABSTRACT 

A  mathematical  model  of  the  Sea  King  Mk.  50  anti-submarine  warfare  (AS W) 
helicopter  and  its  sonar  system  is  presented.  The  model  represents  both  performance 
and  dynamic  flight  behaviour  over  a  range  of  conditions  and  incorporates  the  aero¬ 
dynamics  and  kinematics  of  the  helicopter,  the  control  systems,  pilot  inputs,  the  cable 
and  sonar  dynamics,  plus  wind  and  sea  state  data.  The  aerodynamics  and  kinematics 
model  is  a  three-dimensional  representation  covering  the  operating  flight  envelope, 
where  rotor  aerodynamics  are  based  on  blade-element  and  actuator-disc  theory,  and 
the  control  systems  contain  models  of  both  the  flying  controls  and  automatic  flight 
control  system.  The  sonar  cable  and  trandsucer  mode /  is  formed  by  a  number  of  attached, 
rigid  links,  where  consideration  of  the  forces  on  each  link  enable  the  three-dimensional 
shape  and  motion  of  the  complete  cable  and  transducer  to  be  predicted.  A  description  of 
the  computer  program  for  the  mathematical  model  is  given  and  sample  results  arc 
included. 
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I.  INTRODUCTION 


During  the  past  15  years,  the  Australian  Defence  Science  and  Technology  Organisatio  ■  has 
provided  scientific  assistance  to  the  Royal  Australian  Navy  (RAN)  in  the  area  of  simulation  of 
helicopter  operations.  This  work  began  when  the  RAN  modified  its  Westland  Wessex  helicop¬ 
ters  to  a  new  standard,  which  included  the  fitting  of  a  replacement  sonar  to  improve  anti¬ 
submarine  warfare  performance.1  In  the  type  of  sonar  used,  the  transducer  is  lowered  into  the 
sea  on  a  cable  suspended  from  the  helicopter  when  in  hovering  flight  (known  as  the  cable  hover 
mode  of  operation — Fig.  1)  and  submarine  detection  information  from  the  transducer  is  routed 
via  the  cable  to  processing  equipment  in  the  aircraft.  In  order  to  obtain  optimum  performance, 
the  transducer  must  be  kept  still  and  upright  in  the  water.  To  do  this,  the  plan-position  of  the 
top  of  the  cable  is  altered  through  the  helicopter  motions  and  this  results  in  transducer  move¬ 
ments  being  controlled.  Consequently,  if  plan-position  is  adjusted  to  maintain  cable  angle  at 
the  suspension  point  at  a  predetermined  value  depending  on  the  wind  velocity,  then  the  trans¬ 
ducer  should  remain  still  and  upright  in  the  water  under  all  wind  conditions.  For  effective 
operation,  automatic  control  of  the  helicopter  is  required,  where  cable  deviation  angles  from  the 
vertical  in  pitch  and  roll  are  sensed  at  the  suspension  point  and  used  as  inputs  to  the  aircraft 
flight  control  system. 

A  problem  which  arose  during  the  Wessex  modernization  programme  was  the  instability 
of  the  aircraft  in  the  cable  hover  mode.  To  investigate  this,  a  mathematical  model  of  the  complete 
aircraft,  control  system  and  cable  was  developed1  which  enabled  prediction  of  the  dynamic 
behaviour  of  the  complete  system  in  this  mode.  The  model  was  used  to  solve  the  problem,  which 
stemmed  from  incorrect  gain  settings  in  the  automatic  flight  control  system.  In  addition,  the 
behaviour  of  the  aircraft  part  of  the  model  was  considered  sufficiently  realistic  for  a  range  of 
flight  conditions  that  other  manoeuvres,  such  as  low  rate  dynamic  responses  and  automatic 
transitions  involving  flight  from  cruise  to  hover  and  vice  versa,  could  be  investigated. 

Because  of  the  success  of  this  work,  the  RAN  asked  that  a  model  of  the  Sea  King  Mk.  50 
helicopter  be  developed  when  this  aircraft  and  its  simulator  were  acquired  about  1975.  This 
helicopter  has  a  single,  fully  articulated  main  rotor  of  five  blades  and  a  conventional  pylon- 
mounted  tail  rotor  of  six  blades.  Both  rotors  use  blades  having  a  NACA  0012  aerofoil  section. 
Its  maximum  all-up  weight  is  21,000  Ibf  (mass  =  9545  kg)  and  propulsion  is  by  twin  free- 
turbine  engines  each  having  a  maximum  rating  of  1660  hp  (1239  kW).  The  aircraft  is  fitted  with 
a  Bendix  AN/AQS-13B  sonar  system  for  ASW  work.  To  facilitate  handling  and  allow  transition 
and  ASW  sonar  search  manoeuvres  to  be  performed,  a  Newmark  Mk.  31  automatic  flight 
control  system  is  fitted. 

While  no  inherent  instability  problems  have  occurred  with  Sea  King,  the  model  is  con¬ 
sidered  to  be  a  useful  tool  for  supporting  aircraft  and  simulator  operations,  and  for  conducting 
simulation  studies  of  proposed  helicopter  modernizations.  For  example,  a  proposal  was  made 
during  1978  to  further  update  the  Wessex  aircraft  which  included  the  fitting  of  a  Sea  King  flight 
control  system.  A  feasibility  study2  was  performed  using  the  Wessex  aircraft  model  combined 
with  the  Sea  King  control  system  model  to  check  dynamic  flight  behaviour  and  the  need  for 
alteration  of  system  gains  and  time  constants.  While  the  Sea  King  model  has  been  developed 
from  the  Wessex,  it  has  more  sophistication  in  the  aerodynamics  representation  to  produce  a 
more  accurate  simulation,  and  additionally  a  completely  revised  control  systems  model. 

Previous  work  on  a  similar  system  to  the  one  described  here  was  undertaken  as  long  ago 
as  1960  by  Plessey  (U.K.)  Ltd.,  using  an  analog  computer.  However,  the  models  used  for  both 
aerodynamics  and  control  systems  were  very  simple  and  only  valid  for  small  perturbations  from 
steady  forward  flight.  More  recently,  many  other  models  for  helicopter  simulation  have  evolved 
(e.g.  Wilcock3,  Austin  and  Vann4),  while  some  descriptions  of  flight  control  systems  for  ASW 
helicopters  have  also  appeared  (e.g.  Collomosse,5  Snell ing  and  Cook6).  The  model  described 
here,  which  is  a  development  of  the  Wessex  model  made  by  Packer,1,7  9  is  a  versatile  dynamic 


model  of  the  Sea  King  Mk.  50  helicopter,  its  control  system  and  its  cable,  designed  for  use  on  a 
digital  computer  and  programmed  in  CSMP-IO(ARL)  simulation  language.  While  it  has  not 
yet  been  validated  against  flight  test  data,  a  series  of  flight  trials  has  been  conducted  specifically 
for  this  purpose  using  an  instrumented  Sea  King  helicopter.10  A  report  covering  validation  by 
the  comparison  of  model  and  flight  test  results  using  analog  matching  techniques  will  be  pub¬ 
lished  at  a  later  date.  The  flight  test  data  will  also  be  used  for  the  development  of  procedures 
for  applying  parameter  estimation  methods  to  the  mathematical  modelling  of  helicopters  in 
general. 

2.  MODEL  FORMULATION 

A  block  diagram,  showing  the  relationships  between  the  major  components  of  the  complete 
model  and  specifying  the  inputs  and  outputs  from  each  component,  is  given  in  Figure  2.  Data 
for  the  aerodynamics  were  supplied  by  Westland  Helicopters  Ltd.  and  enabled  a  detailed  aero¬ 
dynamics  and  kinematics  model  to  be  formulated,  based  on  the  physical  theory  established  in 
Reference  11.  Data  and  information  for  the  control  systems  were  obtained  from  a  variety  of 
sources,  including  tests  on  aircraft  equipment,  manuals,  circuit  diagrams  and  data  on  funda¬ 
mental  control  laws  supplied  by  Louis  Newmark  Ltd.  The  data  for  the  cable  and  sonar  transducer 
were  obtained  from  wind-tunnel  tests  on  a  scale  model  of  the  transducer  and  from  information 
supplied  by  Bendix.  This  model  uses  the  theory  developed  in  Reference  8. 

The  aerodynamics  and  kinematics  model  is  a  three-dimensional  dynamic  representation 
where  the  rotor  aerodynamics  are  based  on  blade-element  and  actuator-disc  theory,  rather  than 
detailed  wake  modelling  techniques.  In  this  way  a  realistic,  manoeuvrable  and  versatile  model, 
which  does  not  require  excessive  computing  time  for  its  operation,  is  formulated. 

The  aircraft  control  systems,  which  comprise  a  number  of  interconnected  mechanical, 
hydraulic  and  electronic  subsystems,  are  modelled  in  two  parts— the  flying  controls  and  the 
automatic  flight  control  system  (AFCS).  The  latter  can  itself  be  subdivided  into  an  autostabilizer- 
autcpilot  mode  and  an  ASW  mode.  Both  individual  components  and  overall  systems  are 
modelled  so  that  the  characteristics  of  the  control  systems  are  adequately  represented. 

In  the  cable  and  sonar  model,  the  cable  and  transducer  are  represented  by  a  number  of 
attached  rigid  links,  and  the  effect  of  the  cable  touching  the  funnel  rim  is  included.  The  model 
enables  dynamic  simulation  of  the  system  behaviour  during  sonar  dunking  operations,  but  does 
not  take  acoustic  performance  or  lowering  and  raising  operations  into  account.  Quantities 
obtained  from  the  model  at  the  cable  suspension  point  provide  inputs  to  both  the  aerodynamics 
and  kinematics,  and  control  systems  models. 

The  pilot  model  shown  in  Figure  2  comprises  logic,  switches  and  variable  functions.  The  logic 
and  switches  are  used  to  actuate  the  various  facilities  of  the  control  systems,  with  functions 
being  used  for  trim  settings,  and  pilot’s  stick  and  pedal  movements.  No  attempt  is  made  to  simulate 
pilot  response  as  part  of  the  control  loop. 

In  order  to  model  the  aircraft,  control  systems  and  cable  behaviour  under  operational  con¬ 
ditions,  provision  is  made  to  include  weather  conditions  in  the  simulation.  Both  steady  and 
gusting  winds  in  any  direction  relative  to  the  earth  can  be  specified,  together  with  specifications 
for  the  amplitude  and  frequency  of  surface  waves,  and  the  velocity  vs.  depth  relationships  of 
sea  currents. 


3.  AERODYNAMICS  AND  KINEMATICS  MODEL 

The  model  is  a  dynamic,  three-dimensional  representation  of  body  motion  with  six  degrees 
of  freedom  (see  Fig.  3).  Main  and  tail  rotor  flapping  motions  are  calculated  but  the  rotor  is 
considered  quasi-static,  in  that  the  disc  is  assumed  to  respond  instantaneously  to  the  motion  of 
the  fuselage.  The  inclusion  of  cross-coupling  terms,  such  as  arise  from  angular  rates,  means  that 
the  model  is  not  limited  to  small  perturbation  studies.  It  is  expected  to  be  capable  of  predicting 
low  rate  dynamic  response  manoeuvres,  together  with  performance  characteristics,  within  the 
lower  speed  range  of  the  flight  envelope.  Operation  up  to  an  advance  ratio  (n)  of  0  3  is  of  interest 
but  the  range  of  validity  of  the  model  has  yet  to  be  established  as  it  is  obviously  dependent  on 
the  assumptions  made.  These  assumptions  and  other  features  of  the  aerodynamic  and  kinematics 
model  are  now  discussed  briefly. 
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3.1  Main  Rotor  Aerodynamics 

Using  blade-element  theory,  forces  and  moments  at  the  hub  (origin  of  axes  of  no  feathering— 
ANF  axes)  are  obtained  by  integration  of  blade-element  forces  with  respect  to  radius  and 
azimuth.  The  expression  for  the  normal  air  velocity  component  at  the  blade  element  includes 
helicopter  angular  rate  terms.  The  following  assumptions  are  made: 

(i)  Blades  are  infinitely  stiff  in  torsion  and  bending  and  hence  second  order  flapping  terms 
are  ignored. 

(ii)  Reverse  flow,  blade  stall  and  compressibility  effects  are  absent.  At  the  highest  speed  of 
120  kn  (p  —  0-3).  it  is  likely  that  compressibility  and  stall  effects  become  evident. 

Analysis  of  recent  flight  test  data10  is  expected  to  provide  a  test  of  this  assumption. 

(iii)  An  estimate  of  the  mean  induced  velocity  through  the  rotor  disc  is  obtained  from  the 
well-known  momentum  theory.  Based  on  the  actuator-disc  concept,  the  rotor  is  considered  to 
have  an  infinite  number  of  blades  which  accelerate  air  through  the  disc.  The  rotor  thrust  is  then 
equal  to  the  mass  flow  times  the  total  velocity  increase  in  the  rotor  wake.  According  to  theory, 
the  induced  velocity  at  the  disc  is  half  the  total  increase. 

In  practice  the  mean  induced  flow  is  somewhat  greater  than  the  value  given  by  momentum 
theory  and  an  empirically-based  factor  is  applied.  The  value  of  1-18  used  in  the  Sea  King  flight 
simulator  mathematical  model12  has  also  been  used  in  the  present  work.  In  the  axial  flow  case, 
the  induced  flow  is  assumed  to  be  uniform  over  the  disc.  For  translational  motion,  the  distri¬ 
bution  of  induced  flow,  v.  over  the  disc  is  represented  by  a  Glauert-type  equation  of  the  form 
e  =  I'mean  (I  i  A'(.v  /?)cos  1/1),  where  XIR  and  <f>  are  the  radial  and  azimuthal  co-ordinates  respec¬ 
tively,  and  I'nit-an  is  the  mean  induced  velocity.  This  expression  represents  a  linear  increase  of 
induced  velocity  in  a  streamwise  direction.  No  concensus  exists  as  to  the  appropriate  gradient 
value,  K,  to  be  used,  but  it  is  generally  agreed  to  relate  to  the  wake  sweepback  angle.  x-  In 
Reference  12  a  function  of  x  is  used  which,  in  the  present  work,  is  approximated  by  making  K 
equal  to  x  expressed  in  radians. 

(iv)  Small  angle  approximations  are  used. 

(v)  Blade-element  profile  drag  coefficient  is  independent  of  incidence.  The  complexity 
incurred  in  using  a  quadratic  expression  is  hardly  justified  in  view  of  the  simplifying  assumptions 
made  regarding  the  distribution  of  induced  flow  over  the  disc.  A  constant  value  of  0  012  has 
been  chosen.  This  is  approximately  50",,  greater  than  the  minimum  drag  coefficient  of  a  two- 
dimensional  NACA00I2  aerofoil  section  at  low  Mach  numbers  and  from  a  performance  stand¬ 
point  should  slightly  overestimate  the  required  torque. 

Calculation  of  total  rotor  inflow  and  thrust  coefficient  is  determined  from  an  implicit 
relation  which  is  solved  iteratively.  Within  the  loop  are  allowances  for  blade  tip  losses  and 
ground  effect.  Tail  rotor  thrust  coefficient  is  determined  in  a  similar  manner. 


3.2  Axes  Transformation 

Rotor  forces  and  moment  coefficients  calculated  in  the  ANF  axes  system  (in  this  case  the 
XOZ  plane  contains  the  incident  wind  vector)  undergo  successive  transformation  to  wind  axes, 
shaft  axes  and  body  axes,  so  that  summation  with  fuselage  and  tail  aerodynamic  forces  and 
moments  can  be  made.  Definitions  of  the  axis  systems  used  are  given  in  Reference  7. 


3.3  Aerodynamic  Body  Forces 

Pending  the  acquisition  of  suitable  wind-tunnel  data  relating  to  fuselage  and  empennage 
forces  and  moments,  the  expressions  used  in  the  mathematical  model  of  the  RAN  Sea  King 
simulator12  were  adopted.  While  these  expressions  describe  behaviour  at  high  angler  of  sideslip 
and  attack,  they  are  based  on  limited  wind-tunnel  data  and  their  accuracy  under  these  conditions 
is  uncertain,  but  likely  to  be  adequate  for  low  advance  ratio  flight.  In  order  to  duplicate  the  pitch 
trim  variation  with  speed,  an  empirically  based  weighting  function  modifies  the  pitch  moment 
arising  from  downwash  effects  on  the  fuselage  and  tailplane. 
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3.4  Dynamics  and  Kinematics 

Classical  equations  of  motion,  expressed  in  helicopter  body  axes,  are  written  in  terms  of 
rotor  and  aerodynamic  forces  and  moments  together  with  gravitational  and  cross-coupling 
accelerations.  Using  Euler  angle  notation,  a  matrix  transformation  gives  the  helicopter  velocity 
and  position  in  the  earth  axes  system.  When  combined  with  the  ambient  wind  conditions,  the 
magnitude  and  direction  of  the  incident  wind  vector  is  established  and  hence  the  orientation  of 
the  ANF  axis  system. 

3.5  Control  Inputs 

In  the  case  of  cyclic  pitch  inputs,  allowance  is  made  for  control  phase  angle  and  steady 
state  blade  lag  angle  when  converting  to  actual  blade  angle  changes.  Pitch-flap  and  pitch  lag 
coupling  are  allowed  for  in  determining  the  effective  collective  blade  angle. 

3.6  Rotor  Speed 

As  rotor  speed  is  governed  to  a  nominally  constant  value,  it  was  assumed  at  first  that  speed 
variations  would  be  small  enough  to  be  neglected.  However,  a  preliminary  analysis  of  flight 
test  results10  shows  this  assumption  to  be  invalid  for  manoeuvres  in  which  there  are  large  and 
rapid  changes  in  torque:  e.g.  dynamic  response  tests  with  collective  stick  or  pedal  inputs.  Hence 
an  empirical  rotor  speed  control  model  has  been  developed.  A  block  diagram  of  the  arrangement 
which  comprises  engine  and  rotor  response  characteristics,  fuel  flow  computer,  engine  speed 
governor  and  anticipator  is  shown  in  Figure  4.  Parameters  of  the  model  have  been  adjusted  to 
give  rotor  speed  and  engine  torque  behaviour  similar  to  that  measured  in  flight  tests.  The  model 
is  fully  described  in  Reference  1 1. 


4.  CONTROL  SYSTEMS  MODEL 

For  the  model,  the  control  systems  comprise  the  flying  controls.  AFCS  (both  autostabilizer- 
autopilot  and  anti-submarine  warfare  modes)  and  sensors.  The  overall  configuration  of  these 
components  is  shown  in  Figure  5. 

When  developing  this  part  of  the  model,  it  was  considered  worthwhile  having  two  versions 
available— “full"  and  "simplified".  The  full  version13  15  represents  the  individual  elements 
forming  the  aircraft  system,  while  the  simplified  version16  represents  only  the  overall  control 
laws  and  systems.  Both  are  useful  for  different  types  of  problem  analysis. 


4.1  Flying  Controls 

The  flying  controls  connect  the  pilot  and  AFCS  inputs  to  the  blade  actuating  mechanisms 
as  shown  in  Figure  6.  Both  the  auxiliary  and  primary  jacks  provide  servo  assistance  and  the 
mixing  unit  enables  collective  movements  to  be  cross-fed  onto  cyclic  and  tail  rotor  blade  motions 
to  enhance  the  handling  qualities.  AFCS  signals  are  input  through  Moog  valves  in  the  auxiliary 
servo  unit  which  comprises  four  jacks  (fore-aft.  lateral,  yaw  and  collective),  each  acting  as  a 
limited-authority  series  actuator.  In  addition,  the  unit  provides  supplementary  cyclic  pitch  control 
from  the  beeper  trim  system.  This  enables  fine  adjustment  of  the  cyclic  stick  to  be  made  through 
pilot-operated  trim  switches  and  also  enables  extension  of  authority  for  AFCS  signals  to  occur. 
In  the  collective  and  yaw  channels,  authority  extension  is  provided  by  "open-loop"  spring 
operation  of  the  auxiliary  servos. 

For  the  purpose  of  modelling  both  the  full  and  simplified  control  systems,  it  is  convenient 
to  separate  the  flying  controls  into  cyclic  (fore-aft  and  lateral),  main  rotor  collective,  and  yaw 
(tail  rotor  collective)  channels.  Taking  the  full  model,  each  component  of  the  flying  controls  in 
each  of  these  channels  is  represented  and  an  example  of  a  block  diagram  (the  fore -aft  cyclic 
channel)  is  shown  in  Figure  7.  The  block  diagram  is  formulated  using  the  aircraft  system  and  the 
equations  for  the  mathematical  model  can  be  derived  from  the  diagram.  Modelling  of  the  other 
channels  uses  similar  principles. 
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4.2  Automatic  Flight  Control  System  (AFCS) 

An  AFCS  facilitates  the  handling  of  the  aircraft  and  enables  ASW  manoeuvres,  like  that 
shown  in  Figure  8,  to  be  performed  automatically.  The  control  law's,  which  use  sensor  signals 
as  inputs  and  operate  on  them  to  achieve  the  desired  features,  are  implemented  in  the  electronic 
amplifier  unit  (Fig.  5).  The  electrical  signals  output  by  the  amplifier  unit  operate  the  flying 
controls  via  the  auxiliary  servo  Moog  valves. 

4.2.1  Autostabilizer- Autopilot  Mode 

The  autostabilizer  functions  as  an  attitude  hold  system  by  means  of  which  the  aircraft  is 
stabilized  at  the  pitch,  roll  and  yaw  attitudes  established  through  the  flying  controls.  Autopilot 
facilities  are  heading  hold  and  barometric  height  hold. 

To  illustrate  the  character  of  the  simplified  model,  the  pitch  cyclic  channel  is  shown  in  Figure 
9.  The  attitude  holding  and  stabilizing  characteristics  are  achieved  by  the  use  of  the  pitch  attitude 
signal  supplied  by  the  aerodynamics  and  kinematics  model,  and  its  approximate  derivative.  The 
roll  channel  modelling  is  similar,  while  the  yaw  channel  provides  both  heading  hold  and  rate 
damping.  In  the  collective  channel,  the  height  datum  is  set  when  the  barometric  altitude  hold  is 
enagaged  and  subsequent  deviations  in  height  are  fed  through  the  channel  to  stabilize  the  aircraft 
at  this  datum. 

4.2.2  ASW  Mode 

For  modelling,  it  is  convenient  to  subdivide  this  mode  into  cyclic  channel  and  radio-altitude- 
hold  parts.  The  cyclic  channels  may  be  further  subdivided  into  pitch  and  roll  sections.  In  each 
of  these,  two  control  laws  are  used,  one  for  transition,  doppler  hover  and  air -sea  rescue  manoeu¬ 
vres  (doppler  mode),  and  the  other  for  cable  hover  (cable  mode). 

For  doppler  mode  control  of  the  pitch  channel,  a  law  incorporating  longitudinal  ground- 
speed  error  and  its  integral  is  used.  The  proportional  error  signal  is  the  difference  between  the 
smoothed  doppler  groundspecd  and  its  reference,  which  is  programmed  to  ramp  from  cruising 
speed  to  zero,  or  wee  versa,  over  a  set  time  period  during  transitions.  The  roll  cyclic  channel  is 
similar,  except  that  zero  lateral  groundspecd  reference  is  used. 

The  object  of  the  cable  mode  control  is  to  maintain  the  sonar  transducer  still  and  upright 
in  the  water.  To  do  this  in  calm  conditions,  the  aircraft's  plan  position  is  adjusted  to  minimize 
deviation  in  the  cable  angle  from  vertical  at  the  aircraft.  A  trim  control  is  incorporated  to  offset 
this  angle  for  use  in  windy  conditions,  so  that  the  aircraft  can  maintain  a  position  ahead  of  the 
transducer  to  counteract  the  bowing  effect  of  wind  and  sea  currents  on  the  cable.  The  control 
signal  comprises  the  cable  angle  error  signal,  its  integral,  and  a  damping  signal  derived  from  air¬ 
craft  acceleration.  In  addition  to  modelling  the  control  laws,  other  features  such  as  beeper 
operation  are  incorporated. 

For  radio-altitude-hold  operation,  level  flight  is  maintained  by  setting  the  desired  radio 
height  and  controlling  the  aircraft  to  maintain  this  height  using  radio  altimeter  sensing. 
Operation  during  transition  manoeuvres  is  similar  .except  that  the  demanded  height  varies  linearly 
from  the  aircraft's  cruise  altitude  to  the  hover  height,  or  vice  versa,  over  a  set  time  period. 

5.  SONAR  CABLE  AND  TRANSDUCER  MODEL 

The  representation  of  the  cable  and  transducer  in  the  Sea  King  model  is  basically  the  same 
as  that  described  in  Reference  8  for  the  Wessex  model.  Flowcver.  there  are  some  errors  and 
inadequacies  apparent  in  Reference  8,  and  these  arc  corrected  in  Reference  17.  The  model 
representation  of  the  cable  and  transducer  is  shown  in  Figure  10.  The  cable  is  first  divided  into 
sections  (Fig.  I  On)  whose  lengths,  which  may  be  all  different,  are  specified  by  the  user  of  the 
model,  taking  into  account  the  desired  accuracy  and  execution  speed  for  a  given  simulation. 
The  transducer  is  considered  to  form  the  lower  end  section.  For  each  section,  the  mass  is  assumed 
concentrated  at  the  centre  of  gravity.  Except  for  the  lower  end  section,  this  corresponds  to  the 
section  mid-point.  The  point  masses  arc  then  linked  by  weightless,  rigid  rods,  joined  sequen¬ 
tially  to  one  another  by  frictionless  pivots  (Fig.  I0h).  Gravitational,  fluid  and  tension  forces  are 


considered  in  the  motion  of  each  link  (Fig.  10r)  so  that  a  set  of  simultaneous  differentia!  equations 
can  be  derived.  Small  angle  approximations  are  used  for  the  relative  angular  displacement  of 
joined  links. 

Besides  the  modifications  to  the  Wessex  cable  model  mentioned  above,  a  more  accurate 
representation  of  the  hydrodynamic  characteristics  of  the  transducer,  embodying  data  from 
wind-tunnel  tests  on  a  0-75  scale  m  >del.  has  been  incorporated.  The  tests  provided  aerodynamic 
coefficients  for  drag  and  lift  forces  and  pitching  moment  as  a  function  of  Reynolds  Number  and 
angle  of  incidence.  Examination  of  the  results  showed  that  variation  with  Reynolds  Number 
could  be  ignored  and  that  prov  ided  the  direction  of  fluid  flow  deviates  less  than  30  degrees  from 
lateral  XY  plane,  drag  coefficients  in  the  lateral  and  longitudinal  directions  could  be  represented 
by  constants.  Since  the  model  does  not  include  lowering  or  raising  of  the  transducer,  this  con¬ 
dition  on  relative  fluid  flow  direction  will  be  satisfied  in  most  cases.  In  view  of  the  complexity 
of  representing  pitching  moment  in  the  model  and  its  small  effect  on  the  overall  motion  (i.e. 
generally  small  oscillations  in  transducer  tilt  angle),  it  has  not  been  included. 

Deployment  of  the  sonar  transducer  and  cable  is  made  through  an  opening  in  the  floor 
of  the  aircraft,  know  n  as  the  funnel,  w  ith  the  point  of  suspension  some  distance  above  the  floor. 
The  characteristics  of  this  arrangement  form  part  of  the  model,  including  the  effects  of  cable 
frictional  resistance  on  the  side  of  the  funnel.  Outputs  from  the  cable  model  include  tension 
force  at  the  suspension  point  (supplied  to  the  aerodynamics  and  kinematics)  and  cable  angle 
deviations  from  the  vertical  in  pitch  and  roil  (supplied  to  the  control  systems). 

In  addition  to  the  cable  link  model,  there  have  been  a  number  of  separate  studies  on  various 
aspects  of  the  cable  and  transducer.  In  each  case,  a  number  of  simplifications  have  been  possible 
allowing  a  greater  insight  into  the  effect  of  certain  parameter  changes.  In  one  of  these  studies, 
steady-state  conditions  were  assumed  and  general  solutions  obtained  for  three-dimensional 
motion  in  a  non-uniform  flowlietd.18-  19  These  solutions  are  applicable  to  the  calculation  of 
offset  cable  angle  corrections  for  use  in  cable  hover  manoeuvres  in  windy  conditions.  In  another 
study  on  the  problem  of  increase  in  amplitude  of  oscillation  of  the  transducer  when  winched 
from  the  sea,  the  cable  was  represented  as  rigid  and  weightless,  with  fluid  forces  neglected.20 

6.  OTHER  FEATURES 

6.1  Computer  Program 

The  computer  program  for  the  complete  model  is  written  in  CSMP-IO(ARL).21  which  is 
a  block-oriented  simulation  language.  The  language  comprises  two  parts,  a  modelling  program 
(BOMMP)  and  an  output  program  (TRANS).  The  model  is  expressed  in  coded  form  with  the 
aid  of  block  diagrams  comprising  a  number  of  linked  modules,  called  blocks,  each  one  repre¬ 
senting  a  particular  function  or  operation:  eg.  integration  and  summation.  The  language 
incorporates  “user-defined"  blocks,  written  as  FORTRAN  subroutines,  which  enable  complex 
algebraic  expressions  to  be  handled  conveniently.  A  large  number  of  outputs  may  be  defined 
within  these  subroutines  by  using  “user-output”  blocks. 

In  the  coding,  the  model  is  represented  by  three  types  of  statements— configuration,  para¬ 
meter  and  function.  The  configuration  statements  describe  the  blocks  used  and  specify  the  way 
in  which  they  are  linked  together.  The  parameter  statements  specify  numerical  values  of  para¬ 
meters  associated  with  the  configuration  statements,  such  as  integrator  initial  conditions,  while 
function  statements  specify  the  co-ordinate  pairs  used  to  generate  a  function.  The  output  program 
part  of  the  language  is  capable  of  producing  graphical  and  tabular  results  in  a  variety  of  formats. 

6.2  Sample  Results 

To  demonstrate  some  capabilities  of  the  model,  sample  results  from  three  different  tests  arc 
presented.  In  the  first  test,  the  model  is  programmed  to  perform  the  following  manoeuvre: 

(i)  Initial  condition  -  steady  flight  at  a  forward  velocity  of  90  kn,  altitude  200  ft  (61  m), 
with  autostabili/cr.  neading  hold  and  radio  altitude  hold  engaged. 

(ii)  Transition  to  hover  at  40  ft  <12-2  m)  altitude,  begun  at  5  s. 


(iii)  Dunking  of  sonar  transducer,  begun  at  100  s. 

(iv)  Run  terminated  at  200  s. 


Such  a  manoeuvre  illustrates  the  performance  of  many  parts  of  the  model  and  time  histories  for 
some  selected  variables  are  shown  in  Figure  1 1.  The  main  features  are  summarised  below: 

(i)  Overall,  the  aircraft  performs  the  level  flight,  transition  down,  doppler  and  cable  hover 
manoeuvres  smoothly  and  in  a  reasonable  manner  in  accordance  with  specifications. 

(ii)  Automatic  cyclic  stick  and  pedal  positioning  occur  during  the  run -  these  are  examples 
of  AFCS  authority  extension  via  the  beeper  and  open-loop  spring  systems  respectively. 

(iii)  The  results  shown  are  obtained  directly  from  the  output  program  part  of  the  simulation 
language  (note  that  altitude  is  measured  positive  downwards  because  of  the  (standard) 
conventions  used  for  body  axes  in  the  model). 

The  second  test  (Fig.  12)  illustrates  typical  level  flight  trim  curve  information;  i.e.  the  vari¬ 
ations  of  pitch  and  roll  attitude  angles,  blade  pitch  angles  and  torque  with  airspeed  up  to  220 
ft;s  (130  kn).  For  the  results  shown,  the  aircraft  weight  was  19,000  Ibf  (mass  =  8625  kg)  and  a 
neutral  centre  of  gravity  position  was  adopted.  The  procedure  used  for  setting  up  the  model  in 
trimmed  level  flight  is  similar  to  that  described  in  Reference  17  for  the  Wessex  model. 

The  curves  for  all  variables  have  characteristics  which  are  typical  of  those  for  helicopter 
level  flight  trim  behaviour.  The  curses  of  pitch  attitude  angle  and  fore-aft  cyclic  blade  pitch 
angle  both  exhibit  a  nose-up  hump  in  the  region  20  to  60  ft  s  (12  to  36  kn).  This  occurs  because 
of  the  download  on  the  tailplane  arising  from  sweepback  of  the  main  rotor  wake.  In  the  model 
the  tailplane  force  coefficient  vs.  tailplane  incidence  angle  is  represented  by  the  non-linear 
relationship  specified  in  Reference  12.  In  addition,  the  curses  for  main  rotor  collective  blade 
angle  and  torque  both  have  a  dish  shape  which  is  characteristic  for  helicopters. 

The  third  test  (Fig.  13)  shows  time  histories  for  selected  variables  during  a  typical  dynamic 
response  test,  in  this  case  a  forward  longitudinal  cyclic  step  input  of  0  01  rad  applied  at  an  air¬ 
speed  of  90  kn.  For  comparison  purposes,  response  characteristics  of  the  model  in  both  un- 
stabilized  and  stabilized  (pitch  and  roll  autostabilizers,  height  and  heading  hold  engaged)  con¬ 
figurations  are  shown.  The  stabilizing  action  of  the  attitude  hold  system  is  well  illustrated,  the 
aircraft  adopting  a  steady  pitch  attitude  corresponding  to  the  new  stick  position,  with  the  pitch 
rate  signal  decaying  after  the  initial  action. 


7.  CONCLUDING  REMARKS 

An  aircraft-control  systems-cable  model  of  the  type  described  here  is  suitable  for  investi¬ 
gating  problems  which  may  occur  in  the  use,  modernization  and  performance  assessment  of  the 
Sea  King  helicopter.  Although  flight  trials  have  been  carried  out  to  obtain  a  data  bank  for 
validation  of  the  model,  no  comparisons  of  the  model  w  ith  flight  test  results  have  yet  been  made. 
However,  in  results  such  as  the  samples  shown,  which  illustrate  some  capabilities  of  the  model, 
it  is  seen  to  behave  in  a  reasonable  manner. 

When  considering  where  the  main  model  inadequacies  may  lie.  it  seems  likely  that  the 
assumptions  used  in  modelling  the  rotor  are  the  most  limiting.  This  arises  partly  because  a 
compromise  has  to  be  struck  between  the  amount  of  detail  involved  and  the  amount  of  compu¬ 
tational  time  required  to  solve  the  relevant  equations.  Although  the  mathematical  model 
presented  here  does  not  have  to  run  in  real  time,  as  does  a  flight  simulator  model,  computer  time 
is  still  an  important  factor.  While  the  blade-elemcnt-actuator-disc  representation  is  reasonably 
frugal  in  this  regard  compared  with  detailed  wake  modelling  representations,  the  latter  may  give 
more  accurate  results.  However,  such  methods  are  currently  impracticable  if  manoeuvres  like 
transitions,  where  conditions  arc  continuously  changing  over  reasonably  long  time  periods,  are 
to  be  simulated. 

The  simulation  language  CSMP-IO(ARL)  has  proved  a  viable  means  for  programming, 
the  language  output  program  providing  a  particularly  convenient  means  for  presenting  results. 
While  this  model  and  its  computer  solution  is  not  well  suited  to  real  time  operation,  it  enables 
considerable  detail  to  be  incorporated  and  provides  good  flexibility  for  investigating  a  wide  range 
of  helicopter-control  system-cable  problems. 
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FIG.  1  HELICOPTER  IN  CABLE  HOVER 
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FIG.  2  MODEL  BLOCK  DIAGRAM 
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FIG.  7  FORE-AFT  FLYING  CONTROLS  MODEL 
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FIG.  10  MODEL  REPRESENTATION  OF  THE  CABLE  AND  TRANSDUCER 
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FIG.  1 1  TIME  HISTORIES  FOR  SELECTED  VARIABLES  (TEST  1) 
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FIG.  12  LEVEL  FLIGHT  TRIM  CURVES  FOR  SELECTED  VARIABLES  (TEST  2) 
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FIG.  13  DYNAMIC  RESPONSE  FOR  STEP  CYCLIC  PITCH  INPUT  (TEST  3) 
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